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ABSTRACT 



The production of nitric oxide by soft solar X-rays in 
the E-region of the ionosphere is investigated. An 

empirical expression for the variation in X-ray flux as a 
function of F10.7 is determined. This expression is 
incorporated into a one dimensional diffusive photochemical 
model to compute nitric oxide densities. The results of 
these calculations are compared with NO observations from 
the Solar Mesosphere Explorer satellite. Variations of 
X-ray flux by a factor of 30 over the solar cycle can 
explain the observed variation in nitric oxide densities. 
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INTRODUCTION 



I . 



This thesis investigates the production of nitric oxide 
(NO) by solar soft X-rays in the E-region of the ionosphere 
over the solar cycle. The variation between extreme 
ultraviolet (EUV) flux and X-ray flux as a function of solar 
cycle is investigated as the principal factor influencing NO 
production in the E-region. This research was prompted by 
the hypothesis of Barth et al. (1988 ) that nitric oxide 
variations are produced by variations in the solar soft 
X-ray flux between 20 and 100 A. The main topic will be an 
analysis of the production of nitric oxide at 110 km as 
measured by the Solar Mesosphere Explorer (SME) satellite. 
These measurements are compared to the results of a 
one-dimensional diffusive photochemical model developed by 
Cleary (1985) at the University of Colorado and modified by 
Siskind (1988) to allow for large variations in the solar 
X-ray flux between 18.62 and 50 A. 



A. THESIS GOALS 

This thesis is a theoretical and observational study of 
the production of nitric oxide in the E-region of the 
ionosphere based on variations between EUV and X-ray flux 
over solar cycle. An empirical expression relating X-ray 
and EUV flux is sought in order to better model the 
production and loss of NO in the lower thermosphere. 
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B. THESIS OUTLINE 



The thesis is divided into seven chapters. Chapter II 
presents background information about the composition and 
structure of the ionosphere, and describes some techniques 
for measuring the density profile of various constituents. 
Atmospheric chemistry is reviewed in Chapter III, with the 
emphasis on photochemistry. The continuity equation and 
physical transport mechanisms are discussed, as well as the 
method used by the model in solving the continuity equation. 
The chemistry of nitric oxide is focused on in section C, 
and the concept of X-ray scaling is introduced. Chapter IV 
discusses photochemical modeling and describes the basic 
structure and operation of the photochemical model used in 
this thesis. Chapter V proposes an empirical relationship 
between the X-ray scaling factor (FRACXH) and 10.7 cm flux 
over solar cycle, with results of comparisons between model 
data and SME satellite measurements presented in Chapter VI. 
Chapter VII summarizes the thesis and presents suggestions 
for further research. 



2 



II . 



BACKGROUND 



Nitric oxide makes up only 0.5 parts per billion of the 
Earth's total atmosphere (NOAA, 1976), yet its role as a 
catalyst in the destruction of stratospheric ozone makes it 
one of the most important trace constituents. The 
ionization of NO by geocoronal Lyman-a is largely 
responsible for maintaining the nighttime E-region 
ionization at about 100 km. Additionally, because nitric 
oxide is very reactive, it plays an important role in the 
daytime ionospheric character. Because this thesis concen- 
trates on ionospheric NO production, a brief description of 
the ionosphere is presented. 

A. THE IONOSPHERE 

Above about 60 km, there are sufficient ions and 
electrons to create a region which affects the propagation 
of radiowaves. This region is called the ionosphere . The 
ionosphere is neutral overall. It is composed mostly of 
neutral molecular and atomic gasses which have concen- 
trations that are several orders of magnitude greater than 
the ions and electrons. The variation of electron density 
with altitude leads to three natural subdivisions within the 
ionosphere. These three regions are called the D, E, and F 
regions. The F-region is further divided into two 
subregions called the Fl- and F2-regions. This character- 
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Altitude (km) 



istic layering" of the ionosphere is due to the variation 
of atmospheric composition with altitude and the wide range 
of ionizing wavelengths associated with each neutral 
species. Figure 2.1 shows the different regions of the 
ionosphere . 
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1. The D-Region 

The D-region is the lowest lying ionospheric region 
and is therefore produced by ionizing radiation at those 
frequencies which penetrate deepest into the atmosphere. It 
extends from about 50 to 90 km in altitude. The major 
source of ionization in the D-region is photoionization of 
NO by the Lyman-a 1215.7 A line. 

2. The E-Region 

The E-region lies between about 90 and 130 km, and 
has a daytime electron density about two orders of magnitude 
greater than that of the D-region. Ionization in this 
region is primarily due to photoionization of molecular 
oxygen and nitrogen by EUV and soft X-rays. 

3. The F-Region 

In the F-region, above 130 km, 0, C> 2 , and all 
photoionize, with the primary source above 150 km being 
atomic oxygen. Within the region, the electron density is 
fairly uniform until about 160 km, where it then rises to a 
daytime maximum of approximately 10 /cm at 220-450 km. 
Because of this characteristic profile, the F-region is 
further divided into the Fl- and F2-regions lying at 130-160 
km and 160-600 km respectively. In the Fl-region, ion 
chemistry dominates the characterization of the ionosphere 
whereas in the F2-region physical transport (i.e., dif- 
fusion) is the most important process. 
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B. THE THERMOSPHERE 



The Earth's neutral atmosphere is divided into regions 
primarily on the basis of temperature gradients. Typically, 
at middle latitudes, the temperature falls with increasing 
altitude until reaching a local minimum of around 200 K at 
about 10 km. Here the atmosphere undergoes a temperature 
inversion where the temperature rises with increasing 
altitude until a local maximum of 290 K is reached at about 
50 km. The atmospheric temperature then again falls with 
increasing altitude until an absolute minimum of about 160 K 
is reached at 90 km. This minimum is termed the mesopause 
and marks the beginning of the thermosphere . In the thermo- 
sphere, temperature again increases with altitude, 
approaching a maximum value in the region of 1000-2000 K at 
a height of several hundred kilometers. Above the 
thermosphere, the atmosphere is essentially isothermal. 

C. REMOTE SENSING TECHNIQUES 

Because the ionosphere is a plasma, it will affect the 
propagation of radio waves. It is for this reason that 
radio sounding methods were the first experimental 
techniques used for studying the ionosphere (see e.g., 
Bauer, 1973). Other techniques can be used to obtain 
experimental information about the chemistry and chemical 
processes of the atmosphere. Among the most important of 
these are ground-based observations of airglow emissions and 
the use of satellite-borne photometers and spectrometers. 
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1. Remote Sensing of Chemical Processes 

a. Ground-based Airglow Observations 

The term airglow is used to describe non-auroral 
light emissions from the upper atmosphere. These emissions 
are extremely weak and therefore require very sensitive 
equipment to be detected. A common instrument used for 
airglow measurements is the scanning monochromator which 
detects light in a very narrow bandwidth. A major 
difficulty encountered in airglow detection is the isolation 
of the emission from the strong background of scattered 
sunlight, (see e.g. f McEwan and Phillips, 1975) 

b. Satellite Observations 

Satellite airglow observations offer an 
advantage over ground observations in that there is almost 
no scattered sunlight to obscure the actual airglow 
emissions. An additional advantage lies in the geometry of 
satellite limb viewing (see e.g., Liou, 1980). This makes 
possible the calculation of the column density of an 
emitting species. Through an "onion-skin" inversion method 
the number density of the emitting species can then be 
calculated (see e.g., Cleary, 1985; Hickman et al., 1977). 

2. Radio Sounding of the Ionosphere 

Electron density profiles of the ionosphere can be 
measured using a variable frequency radar which measures the 
time delay between the transmission of an RF pulse and its 
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echo from a reflecting level within the ionosphere (see 
e.g., Bauer, 1973). Such a device is called an ionosonde . 



The requirement for reflection of an electromagnetic 
wave from an ionized medium is that the magnitude of the 
complex propagation vector, K, be imaginary. This require- 
ment is met when 

[ 6 ] - e - ^ > 1 . ( 2 . 1 ) 

e 0 mw 



The value of co for which equation (2.1) is identically one 
is Wp, the plasma frequency. The index of refraction, n, in 
the ionized medium is given by 



2 >. 



n = 






w 



1 - 



80 . 6 [ e ] 



( 2 . 2 ) 



The requirement for reflection of an electromagnetic wave at 
normal incidence is n = 0. Thus the maximum frequency for 
normal reflection is given by 

f 2 = 80 . 6 [ e ] (2.3) 
where f is measured in Hz and [e] is the electron density in 

3 

electrons/m . Thus by varying the frequency it is possible 
to map out the electron density profile. However, 
examination of Figure 2.1 shows that the electron density 
above about 250 km cannot be measured by ground-based 
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ionosondes because the frequencies needed to probe this 
region will already have been reflected at a lower altitude. 

In order to obtain electron density profiles for the 
region above 250 km, ionosondes have been placed on 
satellites. Top-side sounding is a very valuable technique 
for obtaining electron density profiles in the upper 
F2-region. Unfortunately, most of the measurements taken 
from top-side ionosondes do not provide full profiles down 
to the peak altitude so that it is difficult to merge this 
data with that of ground-based ionosondes. 
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Ill . 



CHEMICAL REACTIONS IN THE THERMOSPHERE 



A great variety of atoms and molecules exist in the 
thermosphere. As will be discussed later, much of this 
molecular soup is subject to reaction with solar radiation. 
The resulting photolytically generated species undergo many 
different chemical processes with the surrounding 
atmosphere. This chapter will examine these processes and 
introduce the method used for analyzing the rapidly changing 
composition and structure of the thermosphere. 

A. PRINCIPLES OF PHOTOCHEMISTRY 

Photochemistry is concerned with chemical reactions that 
are initiated by the absorption of solar radiation. The 
absorption process is written as 

A + h\> => A* , (3.1) 

★ 

where A is an excited energy state of molecule A. The 
excited state could be a rotational or vibrational level, an 
electronic state, or it could represent the decomposition of 
the molecule into two or more fragments. This latter 
process is called photodissociation or photolysis if the 
products are all neutral atoms or molecules, and is called 
photoionization if the products include an ion and one or 
more electrons. 
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Once excited molecules or atoms are produced, they 
undergo secondary chemical reactions, forming species that 
would not exist in a normal state of local thermal 

equilibrium (LTE). Therefore these photochemical processes 
are non-LTE mechanisms and must be quantitatively analyzed 
by studying the rate at which reaction occurs (see e.g., 

Chamberlain, 1978 ) . 

1. Photodissociation 

The process by which the absorption of a photon 
leads to disruption of the absorbing molecule is called 
photodissociation or photolysis. The atoms produced in this 
manner are highly reactive. Photolysis processes which take 
place in the photochemical model used for this thesis are: 

N0 + hv=>0 + N( 4 S), (3.2) 

N 2 + hv => N( 4 S) + N( 4 S), (3.3) 

and 

N 2 + hv => N( 2 D) + N( 2 D). (3.4) 

Notice that the last two equations involve the 

photodissociation of molecular nitrogen into two different 
electronic states of elemental nitrogen. The extent to 
which these two reactions occur is given by their branching 
ratios . For the above two reactions, the branching ratio 
for each has been determined to be 0.5, indicating that half 
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and that the 



of the reactions which occur will produce N(^S) 

2 

other half will result in N( D) production (see e.g., 
Heiklen , 1976 ) . 

a. Photodissociation Rates 

There are three factors influencing the rate of 
photodissociation : molecular concentration, photon flux, 

and quantum yield. The most obvious factor must be the 

amount, or concentration, of molecule A. Concentrations of 
a substance per unit volume are denoted by enclosing the 
substance in square brackets. [A] is the concentration of 
substance A per unit volume. 

Another factor influencing the rate is the 
number of photons available to cause photodissociation. The 
flux of solar radiation outside the atmosphere is simply a 
function of the brightness of the Sun and the distance from 
the Sun. However, as the radiation penetrates the Earth's 
atmosphere the flux of photons is attenuated due to 
absorption by atmospheric gasses. At high altitudes the 
atmospheric density is low and therefore little absorption 
takes place. As the altitude decreases the density 

increases, and absorption becomes important. 

The efficiency, or quantum yield of the reaction 

★ 

given by equation (3.1) is defined as the number of A 
molecules produced for each absorbed photon. Quantum yield 
varies with wavelength and is zero longward of the 
ionization wavelength. 
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A common way of expressing the photon flux at a 
particular altitude z is in terms of the exoatmospheric flux 
and the optical depth . Optical depth, x, is defined as 

CO 

t = 2E> a (i)J ( i ]dz ' , (3.5) 

• z 

1 

where a is the absorption cross section for the i*"* 1 

constituent at a specific wavelength, and [ij is the 

concentration of that constituent. The photon flux is then 
given by 

*(z) = *„e -T ’ sec(X) , (3.6) 

where $(z) is the photon flux at altitude z, # is the exo- 
atmospheric photon flux, x is optical depth, and X is the 

solar zenith angle. Note that x*sec(X) is the optical 
pathlength through the atmosphere. 

2. Photoionization 

When the products of equation (3.1) include molec- 
ular or atomic ions and one or more electrons, the process 
is called photoionization. Schematically this is written 

A + h\> =* A + e . (3.7) 

In this process, the energy carried in by the photon is 
converted into ionization energy plus kinetic energy of the 
electron. These energetic electrons are called photo- 
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electrons , (e*), and enter into several different energy 



transfer reactions with neutral gasses. 

The major source of ionizing radiation in the 
thermosphere is solar radiation. The only wavelengths which 
penetrate to the E-layer are extreme ultraviolet ( EUV ) , 
X-rays, and Lyman-a. 

a. Photoionization Rates 

The rate at which photoionization takes place 
for each constituent is described by its ionization rate . 
For a given constituent, the ionization rate is given by 



where o is the ionization cross section and $ is the solar 

z 

flux at altitude z. (Torr et al., 1979) 

B. DIFFUSION AND THE CONTINUITY EQUATION 

It is possible to construct an equation showing the time 
variation of a particular constituent's concentration once 
the rate coefficients for the production and loss mechanisms 
are known. This equation is called the continuity equation 
and contains terms relating the chemical creation and 
destruction of the constituent, as well as terms describing 
the flux due to physical transport processes. There are 
four types of physical transport mechanisms: molecular 

di f fusion , ambipolar diffusion , eddy diffusion , and thermal 
di f fusion. 




(3.8) 
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1. Molecular and Ambipolar Diffusion 

When the concentration gradient of a constituent in 
the atmosphere differs from its equilibrium value, the 
resulting motion is called molecular diffusion. For an 
isothermal atmosphere, this can be written 



*. 

1 




r ati] . [il l 

l 32 + H.J 



(3.9) 



where is the flux, is the molecular diffusion coef- 

t h 

ficient, and is the scale height for the i c constituent. 

Diffusion which results from a concentration 
gradient of a charged constituent in a partially ionized 
gas, or plasma, in which the constituent is influenced by 
the electrostatic forces of the plasma as a whole is called 
ambipolar diffusion . Although electrons in the thermosphere 
have a greater diffusion coefficient than the ions, 
electrostatic attraction between electrons and positive ions 
prevents any large scale separation. Thus the plasma 
diffuses as a whole, with ambipolar diffusion coefficient 
D . (see e.g., McEwan and Phillips, 1975) 

cl 

2. Eddy Diffusion 

Eddy diffusion is the result of turbulent mixing in 
the presence of a composition gradient. Whereas molecular 
diffusion results from the average velocity of the 
constituent being different from that of others, eddy 
diffusion occurs whether or not there is any difference in 
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average molecular velocity between the constituents. Eddy 
diffusion can be expressed similarly to molecular diffusion 
by an equation analogous to (3.8). The coefficient of 
molecular diffusion is replaced by the eddy diffusion 
coefficient K. 

3. Thermal Diffusion 

Thermal diffusion is due to large temperature 
gradients in a gas mixture. It can be properly accounted 
for only in a rather complex manner. However, since the 
thermal diffusion factor is proportional to the 
difference in molecular weights of the constituent gasses, 
thermal diffusion is important in the atmosphere only for 
light gasses such as hydrogen and helium. Because the model 
used in this thesis is primarily concerned with the 
chemistry of nitric oxide, thermal diffusion can be 
neglected . 

4. The Continuity Equation 

Physical transport and chemical processes act on 
each constituent of the atmosphere to produce a variation in 
concentration with time. The equation which governs these 
processes is the mass continuity equation 

^r 1 ♦ - p i - hiii- < 3 - 101 
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Again, [i] is the concentration at time t and 
altitude z, and L^[i] are chemical production and loss 
rates, and is the mean vertical velocity due to transport 
processes. Since [i]*V^ is the vertical flux the second 
term in equation (3.9) may be written 













(3.11) 


For a 


neut 


ral constituent with low 


concentration , 


n may be 


written as 


the sum of molecular and 


eddy 


diffusion, 


, Thus 


= 

i 


= -D. 

l 


ra f i 1 , [i] 3T [in 

taF + T 3 Z + H.J K 


ra [ i ] . [i] 3T 

l 3z T 3 z 


+ [ill 
+ h y 

av 












(3.12) 


where 


the 


middle terms are due to 


the 


dependence 


of scale 


height 


on 


temperature gradient, 


and 


K are the 


molecular 



and eddy diffusion coefficients respectively, T is the 

t h 

absolute temperature, the scale height for the 1 

constituent, and H the scale height for a constituent 
having the mean molecular mass. (Colgrove et al., 1966) 
a. Time Constants 

The continuity equation includes chemical and 
dynamical processes. The dominance of one process over 
another can be determined by examining the relative 

magnitude of the characteristic times, or time constants , of 
the competing processes (Goody and Walker, 1972). If two 
processes are competing and one acts more rapidly than the 
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other, the faster one will be more effective. The dominant 
process will be the one with the shorter time constant (see 
e.g., Wallace and Hobbs, 1977). 

By performing dimensional analysis on the 
continuity equation, the time constants for chemistry, x , 
eddy diffusion, x v , and molecular diffusion, x_., can be 
approximated (Cleary, 1985). The distribution of [ i ] is 
determined principally by chemistry if x~ < x v and x n < x^ , 
and by dynamics if x R < x c and x D < x c . 

A comparison of the lifetimes for molecular 
diffusion, eddy diffusion, and chemical destruction was done 
by Cleary (1985) for the different species used in that 
model. The comparison showed that dynamical processes are 

4 

only important for two species, NO and N( S). For the 
remaining species (N( 2 D), NO + , C> 2 + , N 2 + ' 0+ ' an< ^ 0(^D) + ), 
chemical processes are sufficiently fast that transport may 
be neglected (Cleary, 1985). Ambipolar diffusion for 0 + was 
later added to the model. 

5. Solutions to the Continuity Equation 

For constituents which have sufficiently short 
chemical time constants, the dynamical processes of 
diffusion may be neglected. Thus the flux term in equation 
(3.10) vanishes and the continuity equation reduces to 

= P - L ( i ] (3.10a) 
for N( 2 D), NO + , 0 2 + , N 2 + , 0 + , and 0( 2 D) + . 
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For those constituents which have short dynamical 

4 

time constants (i.e., NO and N( S)), the flux term in the 
continuity equation must be retained. When the expression 
for $ given by equation (3.12) is substituted into equation 
(3.10) and differentiated, a parabolic differential equation 
results. This equation can be written 



am 

at 



iim 
1 a 




3(j] 
i a 



+ 



C.[i] 



+ E i- 



(3.13) 



The coefficients A^ , , C^ , and in equation 
(3.13) are given in Table I. Equations (3.10a) and (3.13) 
cannot be solved analytically and must be solved by 
numerical methods. In the model used for this thesis, the 
method is the Crank-Nicol son finite differencing technique. 
(Cleary, 1985 App. A) 



C. CHEMICAL REACTIONS IN THE LOWER THERMOSPHERE 

Many types of chemical reactions are continually taking 
place in the atmosphere. Products of these reactions include 
photolytically generated atoms, radicals, excited molecules 
and ions. As was seen in previous sections, photoionization 
and photodissociation are the primary photochemical 
reactions. Reactions involving the products of these two 
types of reactions with themselves and with the stable 
atmospheric gasses determine in great part the 
characteristics of the atmosphere. 
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TABLE I 

COEFFICIENTS TO THE PARABOLIC DIFFERENTIAL EQUATION 




The appendix summaries the types of chemical processes 
which are important in the execution of the one dimensional 
model used in this thesis for the calculation of nitric 
oxide in the thermosphere. Of course, many other reactions 
also occur, but it is the chemistry of nitric oxide which is 
germane to this investigation. 
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1. Nitric Oxide Chemistry 

Direct sources of nitric oxide in the E-region are 

4 

the reactions of atomic nitrogen (ground state N( S) and 

2 2 

excited states N( D) and N( P)) with molecular oxygen. 

These three reactions (R28, R29, and R39 of Appendix) do not 

equally contribute to the production of NO. In steady 

state, the NO production rate of each reaction is found by 

multiplying the rate coefficient by the product of the 

concentrations of the reactants. Figure 3.1 shows some 

typical daytime nitrogen density profiles. Figure 3.2 shows 

the corresponding NO production rates from the three 

different atomic nitrogen states, as well as the rate 

coefficient (k) for each reaction. It can be seen from 

Figure 3.2 that the principal source of NO at 110 km is 

2 

N ( D) via the reaction 

N( 2 D) + 0 2 «► NO + 0 (3.14) 

4 

The small contribution from N( S) is due to the temperature 

dependence of the rate coefficient, and the small 

2 

contribution from N( P) is because of its low concentration. 

2 

Sources of N( D) are from the dissociative 
recombination of and NO + (R19 and R27), the dissociative 

excitation by photoelectrons of N 2 (R12), the photo- 
dissociation of N 2 (R15), and the ion-atom interchange 
between and 0 (R17). In a method analogous to that used 

for NO production, it can be shown that the two principal 
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Altitude (km) 




2 

sources of N( D) are via R27 and R12. These equations are 
stated here for ease of reference 

N0 + + e * 0 + N( 2 D) R27 

N 2 + + e* * N( 2 D) + N( 2 D) + e R12 

The two loss processes for nitric oxide in the E-region are 

4 

via the bimolecular reaction with N( S), (R33), and the 



charge exchange with 0 2 + , (R22) 




N( 4 S) + no =♦ n 2 + 0 


R33 


0 2 + + NO =► NO + + 0 2 


R22 
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Figure 3.2. NO Production Rates 
2 

Thus, while N( D) is a source of nitric oxide at 100 km, 

4 

N( S) is both a source and sink. An analysis of these 
production and loss processes as they pertain to the 
photochemical model is presented in Chapter V. 
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IV. THE PHOTOCHEMICAL MODEL 



A. PHOTOCHEMICAL MODELING 

Photochemical modeling is a useful tool for gaining an 
understanding of the minor constituents in the thermosphere. 
By comparing a theoretical model with observations of some 
constituents for which data are available, one can determine 
the validity of the model and better understand physical and 
chemical mechanisms occurring in the thermosphere. 

"Photochemical models are sensitive to various input 
parameters such as solar fluxes, photoabsorption 
cross-sections, chemical reaction schemes and reaction rate 
constants, and transport parameters (e.g., eddy diffusion 
coefficient) as well as to the method of numerical analysis 
and boundary conditions." (Shimazaki, 1985). Additionally, 
models may be time dependent or time independent. Time 
independent models are usually used to arrive at a 
steady-state condition. A steady-state model does not 
represent the actual state for the constituents whose 
concentrations change as a function of time. Causes for 
these changes would be solar flux, geomagnetic Ap index, and 
any temperature dependent reaction rate since the temp- 
erature is also time dependent. Time dependent models 
attempt to account for these temporal variations. 
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B. THE TIME DEPENDENT PHOTOCHEMICAL MODEL 



1. Numerical Technique 

In Chapter III it was shown that the three phenomena 

which must be understood in order to model an atmospheric 

constituent properly are chemistry, eddy diffusion, and 

molecular diffusion. Additionally, it was shown that the 

continuity equation reduces to equation (3.10a) for species 

which do not require transport processes to be considered. 
4 

For NO and N( S) transport processes must be included and 
the continuity equation becomes a parabolic differential 
equation (equation (3.13)). 

Equation (3.10a) can be solved by the following 
semi-implicit finite difference equation 



(U ) t+1 - m t > 

At 



- p t - L t m t + i ■ 



(4.1) 



which is solved for [i]^ + ^, where [i]^. + ^ -*- s the density of 
t h 

the i constituent at time interval t+1 and is a function 
of t i ] t . 

Equation (3.13) is first order in the time 
derivative and second order in the spatial derivative and 
therefore requires three boundary conditions: one in time 

and two in space. Cleary (1985) has shown and Siskind 
(1988) has verified that the spatial boundary conditions are 
satisfied by specifying conditions at the upper and lower 
altitude limits of the model. Several different spatial 
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conditions 



can be used to implement these boundary 

conditions. The ones used by the model require that there 

4 

be zero flux of both NO and N( S) at the upper boundary, 

and at the lower boundary the condition for NO is that there 

4 

is a local minimum in the density profile, while for N( S) 
the requirement is that it be in a state of photochemical 
equilibrium . 

When equation (3.13) is cast into its finite 
difference form using the Crank-Nicolson method, a series of 
coupled algebraic equations are produced which can be 
expressed as a tri-diagonal matrix equation. The density 
profile for each time increment is obtained using the Thomas 
algorithm of matrix inversion. (Von Rosenberg, 1969; Cleary, 
1985) 

2. Inputs to the Model 

The Cleary time dependent photochemical model is run 
in three distinct steps. The first step uses the MSIS-83 
program (Hedin, 1983) to calculate the neutral atmosphere. 
The second step uses a modification of the NRL photoelectron 
code PEGFAC (NRL Memorandum Report 5004, 1982) to calculate 

photoelectron impact rates. Siskind (1988) modified the 
PEGFAC code to extend the energy grid so as to include soft 
X-rays up to 666 eV (18.62 A). The existing energy grid was 
also updated to reflect more recent data. The third step is 
the actual one dimensional model. The model calculates the 
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time development of each species which enters into the 
chemistry of nitric oxide. 

a. Neutral Atmosphere 

The first step of the photochemical model uses 
the MSIS-83 (Hedin, 1983) code to calculate the neutral 
atmosphere for each model "hour" of every day for which the 
model is run. Results are stored as 183 element column 
vectors. Each "hour" five of these vectors are generated. 
The first contains information on the geophysical parameters 
used. These are (1) the exoatmosphe ric temperature, (2) 
solar zenith angle, (3) the F10.7 value for the previous 
day, (4) the 90 day average F10.7 value, (5) geomagnetic Ap 
index, (6) latitude, (7) longitude, and (8) the X-ray 
scaling factor, FRACXR. The second vector contains the 
temperature profile for altitudes from 86 to 450 km in 2 km 
increments (hence 183 elements). The third fourth and fifth 
column vectors contain the density profiles of 0, C> 2 , and 

, respectively. Thus each "hour" is a 183 x 5 array. 
Each day produces 24 such arrays. 

b. Photoelectron Impact Rates 

The calculation of photoelectron impact rates 
was done using the computer program PEGFAC ( NRL Memorandum 
Report 5004, 1982). Recently, Siskind (1988) revised the 

PEGFAC code to account more properly for the higher energy 
photoelectrons produced by soft X-rays. He expanded the 
original 37 EUV intervals from 50-1050 A (Torr et al., 1979) 
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to include an additional 20 intervals which cover the soft 
X-ray flux from 18-50 A. The fluxes used were measured in 
one of Hinteregger ' s rockets and are available from the 
National Space Science Center in a file labelled SC#21REFW 
( Siskind, 1988 ) . 

EUV fluxes are calculated by performing a linear 

interpolation between the appropriate solar flux periods 

reported by Torr et al. (1979) based on the value of FlO.7 
-22 -2 -1 

(xlO W m Hz ). Thus the EUV flux is scaled according 
to Fl 0 . 7 . Since only one measurement of X-ray flux is 
available, interpolation is not possible. However, in order 
to explain the observed solar cycle variations of NO, Barth 
et al. (1988) suggested that the soft X-ray flux exhibits a 
larger variation with solar activity than EUV. For this 
reason, Siskind added a multiplicative X-ray flux scaling 
factor, FRACXR , to the PEGFAC code. The X-ray fluxes used 
by the model are scaled by the multiplicative factor FRACXR . 
It is the goal of this thesis to determine the relationship 
between FRACXR and FlO.7. 

PEGFAC uses the flux data to calculate the 
photoionization rates of N£ , O 2 , and 0, and the photo- 
dissociation rates of N 2 for each model "hour" of every day 
for which the model is run. As with MSIS-83, PEGFAC stores 
the results as 183 element column vectors. For each "hour", 
four of these vectors are generated. The first three 
contain the photoionization rates for N 2 , C> 2 , and 0, 
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respectively. The fourth contains the N 2 photodissociation 
rates. Thus each hour is a 183 x 4 array. Each day 
produces 24 such arrays. 

c. Direct Inputs to the Model 

The neutral atmosphere, photoionization rates 
for N 2 , 0 2 , and 0, and the photodissociation rates for N 2 

are calculated exterior to the model. The data files 
created by these two procedures are used as inputs to the 
model program. The density profiles for the minor species 
are set to a generic initial condition, and the model is 
allowed to run in a time independent mode for a sufficient 
length of time to allow this initial minor atmosphere to 
reach the steady state conditions which correspond to the 
starting date of the run. The model "solves" the continuity 
equation for each of the species involved in the nitric 
oxide chemistry every 15 minutes using the numerical methods 
previously described. When the sun is within 15 degrees of 
the horizon the atmosphere is rapidly changing, therefore 
the time increment of the model switches to five minutes. 
As with MSIS-83 and PEGFAC , the model stores the density 
profile for each minor constituent in a 183 element column 
vector . 

3. Outputs From the Model 

The program is written to calculate the density 
profile for the minor species at 15 minute intervals, or at 
five minute intervals if the sun is within 15 degrees of the 
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horizon. For the purpose of this thesis, the output has 
been structured to extract the nitric oxide density at 110 
km for each day at 1500 local time. This is because the SME 
satellite data was taken at an altitude of 110 km at 1500 
local time each day. The output is stored in a vector 
containing N elements, where N is the number of days for 
which the model was run. 
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V. ANALYSIS 



Chapters III and IV introduced the theory and modeling 
technique that is used in this thesis. This chapter will 
present the experimental approach used to determine the 
relationship between FRACXR and F10.7. 

A. SME OBSERVATIONS 

Data from the Solar Mesosphere Explorer satellite was 
provided by Dr. Charles A. Barth of the Laboratory for 
Atmospheric and Space Physics, University of Colorado 
(telephone conversation 17 November, 1988). Equatorial 
latitude nitric oxide data was used for this thesis. This 
data is the result of averaging measurements from two 
satellite orbits between 20° north latitude and 20° south 
latitude. The first orbit took measurements from 20°-60° 
east longitude, while the second were between 105°-147° west 
longitude. All measurements were taken at an altitude of 
110 km at 1500 local time. The 110 km nitric oxide density 
for the period January 6, 1982 to April 19, 1985 is shown in 
Figure 5.1. 

1. Analysis of Original Model 

Figure 5.2 shows a comparison between SME data and 
model calculations prior to the addition of Siskind's FRACXR 
parameter. It clearly shows that the model was inadequate 
in predicting both the absolute magnitude and variations in 
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Figure 5.1. Nitric oxide density at 110 km. 

SME measurements from January 6,1982 to 
April 19, 1985. 



the nitric oxide density which were observed by the SME 

satellite. The reason for this inadequacy can be explained 

by taking a closer look at the relevant photochemistry. 

Originally, an increase in the EUV flux would increase the 

2 

rate at which the production of N( D) occurs (see R12 and 

R27 in Appendix). However, this increase in EUV flux would 

also cause one of the principal loss reactions for nitric 

oxide (R22) to increase its rate due to an increase in the 

2 

ionized molecular oxygen density. The increase in N( D) 
(and hence NO) was thereby partially offset by the increased 
loss of NO via R22. A process was sought which would allow 
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Date 
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Figure 5.2. Comparison between SME observations 
and the photochemical model prior to FRACXR. 



the model to increase the production of N( D) without 

causing a corresponding decrease in NO density via the 

loss reaction R22. This led to the modification of the 

PEGFAC code by Siskind as discussed in Chapter IV. 

The reason behind the incorporation of the 

multiplicative X-ray scaling factor is now clear. Soft 

X-rays which penetrate to the E-region photoionize the major 

constituents and produce energetic photoelectrons as 

discussed in Chapter III. The increase in photoelectrons 

2 

drives the production of N( D) via R12 and R17, but does not 
affect R22. In order for this mechanism to work, the 
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variation in solar X-radiation must be greater than the 
variation in solar EUV radiation. (Barth et al., 1988) 

B. DETERMINATION OF FRACXR 

Initial estimates of the magnitude of FRACXR were 
provided by Siskind (telephone conversation November 9, 
1988). He suggested that FRACXR take on values near one for 
conditions of solar minimum and 50 for solar maximum. This 
provided a starting point for the determination of the 
relationship between FRACXR and solar cycle. 

1. Assumptions 

Two major assumptions were made for the purposes of 
simplifying the data analysis and reducing computer time. 
The first assumption is that auroral electrons do not 
influence the nitric oxide chemistry. This assumption is 
based on the fact that all computations and observations 
were at equatorial latitudes. It is justified by the work 
of Siskind (1988), who showed that for equatorial latitudes 
there is little enhancement of NO density due to an auroral 
storm. 

The second assumption was that the computation of 
nitric oxide density would be independent of the longitude 
at which the model was run. This allowed the model to be 
run at only one longitude for the computation of each NO 
density time sequence, instead of averaging two sets of 
measurements as was done for the SME data. To justify this 
assumption, the model was run at both 40 degrees east 
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NO Density (xlO 7 cm ' 



longitude and 126 degrees west longitude for the same five 
day period in 1982. Figure 5.3 shows the result of this 
procedure. The two cases deviated from each other by less 
than 0.2%. Based on this result, all subsequent model runs 
were done only at 40 degrees east longitude. 




Date 

Figure 5.3. Model [NO] calculations for 
two different longitudes. 

2. Experimental Approach 

The SME data was searched to locate periods where 
the value of F10.7 remained constant for at least three 
days. Four such periods were initially found corresponding 
to F10.7 ( xlO -22 Witf 2 Hz -1 ) values of 70, 90, 128, and 201. 
The model was then run with various values of FRACXR for 
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each of the four periods. Figures 5.4 through 5.7 show the 
results of these runs. This initial data is summarized in 
Table II and shown graphically in Figure 5.8. Uncertainties 
are estimated from the data directly without measuring them 
explicitly since only one set of calculations is possible 
each time the model is run (see e.g., Bevington , 1969 ) . 





SUMMARY 


TABLE II 

OF INITIAL FRACXR DATA 




DATE 

START END 


F10.7 

(xlO -22 Wm“ 2 Hz“ 1 ) 


FRACXR 


84300 


84303 


70 


3 ± 2 


83297 


83300 


90 


6 ± 1 


83185 


83188 


128 


22 ± 2 


82169 


82171 


201 


24 ± 1 



It can be seen from Figure 5.8 that the first three 
data points are nearly linear, while the point corresponding 
to an F10.7 value of 201 lies significantly lower than would 
be expected using a linear relationship based on the first 
three points. To resolve this uncertainty, a fifth data 
point was taken for an F10.7 value of 242 . Although this 
value was constant for only two days, the FRACXR value at 
this point would determine if a linear relationship was 
probable. The result of this run is shown in Figure 5.9. 
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Figure 5.4. Model [NO] densities for F10.7 = 70. 
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Figure 5.5. Model [NO] densities for F10.7 = 90. 
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Figure 5.6. Model {NO] densities for F10.7 = 128. 
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Figure 5.7 
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Model [NO] densities for F10.7 = 201. 
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Figure 5.8. Initial FRACXR Data. 
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Figure 5.9. Model [NO] densities for F10.7 = 242 
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3. Fitting the Data 

The addition of the 
is reflected in Table III, 
data. The problem now is to 
of functions) which best fi 
then be used to determine 
basis for input to the model 



fifth data point at F10.7 = 242 
a summary of the final FRACXR 
find a function (or combination 
ts the data. This function can 
the value of FRACXR on a daily 





SUMMARY 


TABLE III 

OF FINAL FRACXR DATA 




DATE 

START END 


F10.7 

, , n -22 ri -2 tI -1 . 
(xlO Wm Hz ) 


FRACXR 


84300 


84303 


70 


3 ± 2 


83297 


83300 


90 


6 ± 1 


83185 


83188 


128 


22 ± 2 


82169 


82171 


201 


24 ± 1 


82342 


82343 


242 


27 ± 1 



Several fitting functions were tested. The criteria 

for determining the best fit was that function which 
2 

minimized x and was monotonically increasing in FRACXR. 

2 

"The definition of X as the ratio of the estimated variance 
2 2 

s to the parent variance a ...makes it a convenient measure 
of the goodness of fit." (Bevington, 1969). A weighted 
least-squares fit to a third degree polynomial is shown in 
Figure 5.10, and a weighted least-squares quadratic fit to 
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the first four points is shown in Figure 5.11. Weighting 
was determined for each data point by 



w . 
1 




/ 



where was taken to be the estimated 
5.12 shows a least-squares logarithmic 
shows a combination of two linear 
Results are summarized in Table IV. 



uncertainty. Figure 
fit, and Figure 5.13 
least-squares fits. 



TABLE IV 

SUMMARY OF FITTING FUNCTIONS 


FUNCTION 




x 2 


THIRD DEGREE POLYNOMIAL 




10.894 


QUADRATIC 




8.392 


LOGARITHMIC 




87.632 


DOUBLE LINEAR FIRST 


LEG 


3.951 


SECOND 


LEG 


0.749 
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Figure 5.10. Weighted Third Degree Polynomial Fit. 




Figure 5.11. Weighted Quadratic Fit. 
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FRACXR 




Figure 5.12. Weighted Logarithmic Fit. 




Figure 5.13. Weighted Double Linear Fit. 
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4. Statement of Empirical Relationship 

Based on the results of the analysis presented in 
the preceding section, the empirical formula for the scaling 
factor FRACXR which relates X-ray flux and EUV flux in the 
E-region of the ionosphere is 

FRACXR = 14.930 + 0.048x(Fl0.7) 
for F10.7 > 130 and 

FRACXR - -24.513 + 0 . 351x( FlO . 7 ) 

for FlO. 7 < 130, where FlO. 7 has units 10 -22 Wm -2 Hz -1 . 
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VI . RESULTS 



In Chapter V, an empirical expression for FRACXR was 
proposed. In this chapter, the expression is incorporated 
into the one dimensional diffusive photochemical model to 
compute nitric oxide for two 60 day time periods. 

A. COMPARISON OF MODEL AND SME DATA 

Two 60 day periods were chosen to further test the 

model. Each period was chosen to correspond to a time where 

the SME data showed distinct periodic variations in nitric 

oxide density. The first 60 day period was from 82213 to 

-22 -2 

82273, and was chosen because the FlO.7 value ( xlO W m 

Hz - 1) ranged from a minimum of 120 to a maximum of 214 and 

thus would correspond to FRACXR values taken from the upper 

curve of Figure 5.13. The result of this 60 day run as well 

as the daily value for FlO.7 is shown in Figure 6.1. Model 

calculations without using X-ray scaling are included for 

comparison. The second 60 day period was from 84210 to 

-22 -2 

84180 , and was chosen because the FlO.7 value (xlO W m 
Hz - ^) ranged from a minimum of 90 to a maximum of 183 and 
thus would correspond to FRACXR values taken from the lower 



curve 


of 


Figure 


5.13. 


The 


result of 


thi s 


60 day run, along 


wi th 


the 


F10 .7 


values 


and 


unsealed 


model 


calculations, is 


shown 


in 


Figure 


6.2. 
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Figure 6.1. Comparison of SME Data to Model 
for Period of Solar Maximum. 



46 



F10.7 (xlO'° Km' 2 Hz' 1 ) 





V 

a 



o 

* 




Date 




Date 



Figure 6.2. Comparison of SME Data to Model 
for Period of Solar Minimum. 
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1. Discussion of Results 

The average daily deviation of the model calculat- 
ions for nitric oxide density from the SME satellite 
measurements was 13.92% for the first 60 day run and 24.27% 
for the second 60 day run. This is not surprising, since 

the first run was based on a FRACXR relationship which has a 

2 

38% probability of being the correct function (x = .75), 

while the second run was based primarily on a FRACXR 

relationship which only has a 4.8% probability of being the 

2 

correct function (x = 3.95). Figures 6.1 and 6.2 show that 
variations of X-ray flux by a factor of about 30 over the 
solar cycle can explain the observed variation in nitric 
oxide density . 

Calculated nitric oxide densities for both 60 day 
runs are plotted against F10.7 to test the correlation 
between the two data sets in Figure 6.3. An inspection of 
the figure shows that correlation does exist and follows the 
same "double linear" type of relation as exists between 
FRACXR and F10.7. A least-squares fit to each branch of 
this scatter plot yields lines with regression coefficients 
of 0.55 and 0.90 for the upper (F10.7 greater than 130) and 
lower (F10.7 less than 130) parts, respectively. This 
strongly suggests that FRACXR is incorrectly calculated for 
values of F10.7 greater than 130. 

An inspection of Figures 6.1 and 6.2 reveals that 
the derived relationship between FRACXR and F10.7 does not 
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Figure 6.3. Correlation of Calculated NO 
Densities and F10.7. 



the model to correctly calculate NO densities for values of 

F10.7 (xlO - ^ W m - ^ Hz - 1) from about 140 to 160. The 

computed NO densities for these values are too high due to 

2 

the artificially high production rate of N( D) from the 
increased X-ray flux. This suggests that a smoother 
transition exists between the lower and upper curves of 
Figure 5.13 than proposed by the "double linear" fit. 
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VII . 



CONCLUSION 



A. SUMMARY 

A theoretical and observational study of the production 
of nitric oxide in the E-region of the ionosphere has been 
presented. The observational basis was the NO measurements 
made by the Solar Mesosphere Explorer satellite. An 
empirical expression for the variation in soft solar x-ray 
flux as a function of the 10.7 cm solar flux was determined. 
This variation was in the form of a multiplicative scaling 
factor, FRACXR, between X-ray and EUV flux given by 

FRACXR = 14.930 + 0 . 048x( F10 . 7 ) 
for F10.7 greater than 130 and 

FRACXR = -24.513 + 0 . 351x( F10 . 7 ) 

for F10.7 less than or equal to 130, where F10.7 has units 
-22 -2 -1 

10 Wm Hz . The 10.7 cm flux was chosen as the 

independent parameter because it is used in the neutral 
atmosphere model as an indicator of solar EUV radiation. 

The expression which calculates FRACXR was incorporated 
into a one dimensional diffusive photochemical model to 
compute nitric oxide densities. Model calculations for NO 
densities showed significant improvement in both magnitude 
and variation. These model calculations were compared to 
observational measurements from the SME satellite. Results 
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of this comparison showed that model calculations using the 
above X-ray scaling factor were in much better agreement 
with the observed SME data than the calculations made 
without X-ray scaling. The conclusion from this invest- 
igation is that variations of soft solar X-ray flux by a 
factor of 30 over the solar cycle can explain the observed 
variation in nitric oxide densities. 

B. SUGGESTIONS FOR FURTHER RESEARCH 

Figure 7.1 shows a scatter plot of observed nitric oxide 
densities versus 10.7 cm solar flux. Barth et al. (1988) 
pointed out that the observed nitric oxide densities show a 
strong correlation with Lyman-a flux over solar cycle. This 
is shown in Figure 7.2. Examination of these two figures 
suggests a stronger correlation of NO density to Lyman-a 
than to F10.7. It is suggested that further investigations 
be made into model calculations of NO density based on an 
X-ray scaling factor which is a function of Lyman-a. 
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Figure 7.1. Correlation of NO and F10.7. 
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Figure 7.2. Correlation of NO and Lyman-a. 
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APPENDIX 



NITRIC OXIDE REACTIONS IN THE LOWER THERMOSPHERE 



This appendix lists all of the reactions used in the 
Cleary one dimensional photochemical model (1985). 
Reactions between neutral species are listed first, followed 
by ion-neutral reactions. Equation numbers shown here 
reflect the scheme used in Cleary (1986), except for R39 



which was added later. 

A. NEUTRAL SPECIES 

1. Photoionization 

0 + hv «*■ 0 + + e Rl 

O 2 + hv * 02 + + e R2 

N 2 + hv =» N 2 + + e R3 

0 + hv =» 0( 2 D) + + e R4 

N( 4 S ) + hv =» N + + e R8 

O 2 + hv **- 0 + 0 + + e R9 

N 2 + hv N( 4 S) + N + + e RIO 

NO + hv * NO + + e Rl3 

2. Photodissociation 

N 2 + hv * N( 4 S) + N( 4 S) Rl 4 

N 2 + hv =* N( 2 D) + N( 2 D) Rl 5 

NO + hv =» 0 + N( 4 S) R16 
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3 . 



Radiative Relaxation 





N( 2 D) =» N( 4 S ) + hv 


R32 


4 . 


Quenching 






N( 2 D) + e =» N( 4 S ) + e 


R30 




N( 2 D) + 0 » N( 4 S) + 0 


R31 


5 . 


Bimolecular Reaction 






N( 4 S ) + 0 2 NO + 0 


R28 




N( 2 D) + 0 2 NO + 0 


R29 




N( 2 P) + 0 2 ■» NO + 0 


R39 




N ( 4 S ) + NO * N 2 + 0 


R33 




n( 2 d) + no => n 2 + 0 


R34 


6 . 


Energy Transfer or Photoelectron Impact 






* + 

0 + e =* 0 + 2e 


R5 




* + 

0 2 + e ^ °2 + 2e 


R6 




N 2 + e* =* N 2 + + 2e 


R7 




N 2 + e* =* N( 4 S) + N( 4 S) + e 


Rll 




N 2 + e* =* N( 2 D) + N( 2 D) + e 


R12 


IONS 




1. 


Charge Transfer 






N 2 + + o 2 * 0 2 + + N 2 


R18 




0 + + 0 2 * 0 2 + + 0 


R21 




0 2 + + NO ** NO + + 0 2 


R22 




o( 2 d) + + N 2 =♦ n 2 + + 0 


R3 5 




o( 2 d) + + 0 2 «► 0 2 + + 0 


R36 
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R17 



2. Ion-Atom interchange 

N 2 + + 0 «► N0 + + N( 2 D) 

0 + + N 2 N0 + + N( 4 S) R20 

0 2 + + N( 4 S) «► N0 + + 0 R23 

N + + 0 2 =» N0 + + 0 R25 

3. Dissociative Recombination 

N 2 + + e => N( 2 D) + N( 2 D) R19 

0 2 + + e =» 0 + 0 R24 

N0 + + e + 0 + N( 4 S) R26 

N0 + + e + 0 + N( 2 D) R27 
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